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ABSTRACT 

Context. Lyman-alpha (Lya) is now widely used to investigate the galaxy formation and evolution in the high redshift universe. . 
However, without a rigorous understanding of the processes which regulate the Lyo- escape fraction, physical interpretations of high-z 
observations remain questionable. 

Aims. To examine six nearby star-forming galaxies to disentangle the role of the dust from other parameters such as gas kinematics, 
geometry and ISM morphology in the obscuration of Lya. Thereby we aim to understand the hya escape physics and infer the 
implications for high-redshift studies. 

/Wefhods. We use HST/ACS to produce continuum-subtracted Lya maps, and ground-based observations (ESO/NTT and NOT) to 
map the Ha emission and the extinction E(B-V) in the gas phase derived from the Balmer decrement HafHp. 

Results. When large outflows are present, the Lya emission appears not to correlate with the dust content, confirming the role of 
the Hi kinematics in the escape of Lya photons. In the case of a dense, static Hi covering, we observe a damped absorption with a 
declining relationship between Lya and E(B-V). 

We found that the Lya escape fraction does not exceed 10% in all our galaxies and is mostly about 3% or below. Finally, because of 
the radiative transfer complexity of the Lya line, star formation rate based on Lya luminosity is underestimated with respect to that 
derived from UV luminosity. Simple reddening correction does not reconcile SFR(Lya) with the total star formation rate. 
Conclusions. The dust is not necessarily the main Lya escape regulatory factor ISM kinematics and geometry may play a more 
significant role. 

The failure of simple dust correction to recover the intrinsic Lya/Ha ratio or the total star formation rate should prompt us to be more 
cautious when interpreting high-z observations and related properties, such as SFRs based on Lya alone. To this end we propose a 
more realistic calibration for SFR(Lya) which accounts for dust attenuation and resonant scattering effects via the Lya escape fraction. 

Key words. Galaxies: starburst - Galaxies: ISM - Ultraviolet: galaxies - ISM: dust, extinction - Galaxies: individual: (Haro 11, 
ESO 338-04, NGC 6090, IRAS 08339-1-6517, SBS 0335-052, Tololo 65) 



1. Introduction 

The Lyman-alpha emission (Lya) has become the most pow- 
erful tracer of star-formation in the high-redshift universe. It be- 
comes the strongest emission line in the optical and near infrared 
(NIR) domain at redshift z > 2.1, and is likely to remain a very 
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competitive tool, even with the advent of extremely large tele- 
scopes (ELTs) and the James Webb Space Telescope (JWST). It 
is widely used as an efficient detection and redshift confirmation 
tool for distant galaxies, to derive star formation rates (SFRs), as 
well as to probe the ionisation state of the intergalactic medium 
(IGM) at the final stage of the reionisation epoch. In this way, the 
last decade has been the high-redshift era, in which the devel- 
opment of new techniques and facilities have allowed the explo- 
ration of the properties of galaxy population, resulting in a major 
improvement in our understanding of the distant Universe. 

The importance of the Lyo' emission line in the cosmological 
context was predicted very early by Partridge & Peebles ( 1967) , 
who suggested that young high-z galaxies, undergoing their first 
star formation events, should be detectable thanks to their strong 
Lya emission. Unfortunately, the first attempts to detect such ob- 
jects stand in contrast with those predictions. Indeed, initial sur- 
veys ( e.g. [Pritchet & Hartwick|[l989 , Djorgovski & Thompson 



|1992[ |de Propris et al.||1993[ l failed to discover the predicted 



space density of Lya emitters. Unsuccessful campaigns and faint 
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Lya fluxes were attributed to the dust attenuation coupled to 
the resonant scattering of the Lya line (see Pritchet 1994' for 
a review). The first break-through came with the observations 
of |Cowie &"H u (1998) and Hu et al.| ( |1998"l l which encouraged 
high-efficiency surveys for high-z LAEs detection. Most no- 
tably, two techniques are now used to detect high-z galaxies. The 
Lyman Break Technique (Steid el et al.||1996) l uses the absorp- 
tion bluewards the Lya absorption edge to detect the so-called 
Lyman Break Galaxies (LBGs). Narrow-band imaging surveys 
that make use of the Lya recombination line produced by the 
reprocessed absorbed radiation to target Lyman-Alpha Emitters 
(LAEs). The success that has been attributed to these techniques 
has given Lya a key-role in the context of understanding the dis- 
tant universe. Besides its use in the identification of galaxies, it 



is used to put constraints on the cosmic reionisation (Malhotra & 
[Rhoads 2004; Kashika wa etaH20 06, Dijkstra et al. 2007|l, and 
to study the clustering properties and morphology of galaxies 
to the highest redshifts ( Hamana et al.|2004; ,Ouchi et al.,2005, 
Murayama et al.|[2"007| l. Finally, it allows the estimation of the 



star formation rates at high redshift ( Kudritzki et al. 2000 Fujita 
[eFal. 2003 ; Pirzkal et al.|2007| l. 

The Lya star formation rate is typically derived using the 
Ha calibration relation (Kennicutt 1998| l and assuming a case 
B recombination theory for the Lya/Ha line ratio. Nevertheless, 
SFRs inferred from the UV continuum are found to be incon- 
sistent with SFR(Lya). It appears that SFR(Lya) is typically 
lower than SFR(UV) by a factor of several (|Ajiki et al ||200J 



ITaniguchi et al.|2005|[T^ken et al.|2007l|Gronwall et al.|2007| l. 
Correction for internal reddening may, in principle, reconcile 
these two indicators. Yet different extinctions experienced by 
the continuum and the emission line may arise due to geomet- 
rical effects ( jCalzetti et"aL||1994[ [Giavalisco et al.||1996| ), and 
radiative transfer effects of the Lya line make this issue far 
from being resolved. This demonstrates that caution should be 
taken when using such a calibration. This need for caution is 
further demonstrated in cosmological studies where only a frac- 
tion of UV-selected galaxies show Lya in emission ( Shapley| 
|et al. 2003). Furthermore some high-z studies have unveiled very 
high rest frame Lya equivalent widths (EWs,jRhoads et al. 2003 



[Finkelstein et al. 2007) and it seems unlikely that such high EWs 
can result from the ionising output of a normal stellar population. 
In the case of inhomogeneous ISM, where the dust is distributed 
in neutral clouds, with an ionised inter-could medium, Neufeld| 
( |1991| l, and later Hansen & Peng Oh| ( |2006[ ) have shown that 
Lya photons could potentially escape easier than continuum ra- 
diation. In this scenario, intrinsic Lya EWs are expected to be 
enhanced, resulting in the large values observed in those studies. 

Many efforts have been devoted to understanding the physi- 
cal processes that govern the fraction of escaping Lya photons, 
regarding its potential cosmological importance. Early observa- 
tions at low redshift (jMeier & Terlevich|1981 Deharveng et al. 



|1985[ [Terlevich et al.|l993 [) have shown Lya to be much weaker 
than predictions or even absent from starburst galaxies. This 
weakening was first attributed to the dust attenuation, confirmed 
by the correlation observed between EWi,ya and the metallic- 
ity (Chariot & Fall 19931. These interpretations notwithstand- 
ing, TGlavalisco^^er^lJi 1996) reached the opposite conclusion, 
finding no clear correlation between EW^yo or Lya/H/? and E(B- 
V). Furthermore, dust correction failed to recover the intrinsic 
Lya/Fly6 ratio predicted by case B recombination theory. Recent 
spectroscopic studies have outhned the Lya "observational puz- 
zle". |Enffi[eFaI](|T994| and |Thuan & Izotov| (1997), with the 
Goddard High Resolution Spectrograph (GHRS) and the Space 
Telescope Imaging Spectrograph (STIS), have found damped 



Lya absorption in IZw 18 and SBS 0335-052, the most metal- 
deficient galaxies known at low-z. In the purely dust-regulated 
model, a prominent Lya emission feature would be expected. 
On the other hand Leq ueux et al.| ( pl^95| have found strong Lya 
emission in a much more metal- and dust-rich starburst Haro 2. 
Further studies have shed light on the mechanisms by which Lya 
photons may escape their host. Kunth et al. (1998 ) observed Lya 
morphologies ranging from emission to absorption in 8 low-z 
starbursts. They found systematic blueshifts between Lya fea- 
ture and Low lonisation State (LIS) metal absorption features in 
the ISM when Lya is seen in emission, indicative of an outflow- 
ing neutral medium. P-cygni profiles, with a redshifted emission 
peak with respect to the systemic velocity, were also found in 

1(|2003i 



these spectra. Furthermore, Mas-Hesse et al. 
hydrodynamic models ( jTenorio-Tagle et aT 



made use of 
1999[) to propose a 



unified physical picture for the variety of observed Lya profiles 
as a function of the evolution of the burst and the viewing ge- 
ometry. They found that the Lya emission visibility and shape 
are mostly driven by the kinematical configuration of the neutral 
gas. Likewise, over the past few years, theoretical works and nu- 
merical simulations have been developed with the same purpose. 
Ahn et alT| (|20()3[) and |Verhamme et al.| P006| l have shown how 
the variety of Lya profiles are mostly driven by the expansion of 
a super-bubble of neutral gas and the properties of the ISM (Hi 
column density and dust content). Hansen & Peng Oh (20061) uti- 
lized the original idea of Neufeld ( 1991| to investigate the effects 
of a multi-phase ISM. More cosmological-oriented simulations 
(Tasitsiomi 2006) have been carried out, although the effects of 
dust remain to be treated. The simulations are currently able to 
deal with arbitrary intrinsic emission characteristics, hydrogen 



density, velocity fields and dust distributions (Verhamme et al. 



2006 1, and to consistently reproduce the Lya profiles observed 
in z ~ 3 LBGs ( [Verhamme et al.|[2007[ [Schaerer & Verhamme] 
2008). 

The complex nature of the Lya escape probability revealed 
by low-z spectroscopic studies, gave a rise to additional issues. 
The resonant scattering phenomenon of Lya line may cause Lya 
photons to travel and to be emitted far from their production sites 
hence to be spatially uncorrected with non-resonant radiation 
(Ha or continuum photons). UV- targeted spectroscopic studies 
may therefore miss a significant fraction of the Lya emission 
if scattered away from the aperture. Ionised holes in the ISM 
may also allow the escape of Lya photons in a spatially lim- 
ited region, and transmission may vary greatly on small scales 
across the starburst region. These consideration are the motiva- 
tion for our Lya imaging survey with the Advanced Camera for 
Survey (ACS) on board the Hubble Space Telescope (HST). We 
have observed a hand-picked sample of six nearby star-forming 
galaxies, with the aim of exploring a large range of relevant pa- 
rameters. Preliminary results have been presented in Kunth et al. 
(|2003| ), and more de tailed studies on ES O 338-04 d Hayes et aL 
|2005[ ) and Haro 11 (Hayes et al. |2007) ). Emission and absorp- 
tion were found on very small scales in central regions of the 
starburst while absorption is observed in front of many of the 
brightest UV sources. 

In this article, besides the HST observations, we make use 
of ground-based observations from the ESO New Technology 
Telescope (NTT) and the Nordic Optical Telescope (NOT) to 
map Ha and Hfi emission and derive the extinction from the 
Balmer decrement (Ha/Hy6). Investigating correlations between 
E(B-V) and Lya emission enables us to disentangle the role 
of the dust from other parameters in the Lya escape mecha- 
nism, and investigate the implications of proper dust correction 
on high-redshift studies. The paper is structured as follows: In 
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Table 1. Targets general pr operties. E(B-V)mw is the G alactic extinction given 



AB system. References: 1: Bergvall & Ostlin (20021, 2: Papaderos et al. (20061, 3: Gonzalez De 
( |200T] | 



by Schlegel et al. 



1998 I and Mr rn agn itudes are in 



gado et al. (1998 1, 4: Izotov et al. 



Target 


Other 


RA(2000) 


Dec(2000) 


E(B - V)mw 


z 


12+ 


Mb 


Ref 


Name 


Name 










log(0/H) 






Haro 1 1 


ESO 350-38 


00:36:52.5 


-33:33:19 


0.049 


0.020598 


7.9 


-20 


1 


SBS 0335-052 


SBS 0335-052E 


03:37:44.0 


-05:02:40 


0.047 


0.013486 


7.3 


-17 


2 


IRAS 08339+6517 


PGC 024283 


08:38:23.2 


+65:07:15 


0.092 


0.019113 


8.7 


-21 


3 


Tololo 65 


ESO 380-27 


12:25:46.9 


-36:14:01 


0.074 


0.009 


7.6 


-15 


4 


NGC 6090 


Mrlc 496 


16:11:40.7 


+52:27:24 


0.020 


0.029304 


8.8 


-21 


3 


ESO 338-04 


Tol 1924-416 


19:27:58.2 


-41:34:32 


0.087 


0.009633 


7.9 


-19 


1 



Target 


Ha 




Hot continuum 








IIy6 continuum 






Camera/Filter 


Exp. 


Camera/Filter 


Exp. 


Camera/Filter 


Exp. 


Camera/Filter 


Exp. 


ESO/NTT observations 


















Haro 11 


EMMI-R/598 


900 


EMMI-R/597 


1200 


SuSI2/549 


2866 


EMMI-R/770 


1800 


SBS 0335-052 


EMMI-R/598 


1800 


EMMI-R/596 


4800 


SuSI2/548 


900 


EMMI-R/771 


4500 


ESO 338-04 


EMMI-R/597 


1800 


SuSI2/778 


3600 


SuSI2/719 


2400 


EMMI-R/771 


4500 


Tololo 65 


SuSI2/709 


1200 


SuSI2/778 


1200 


SuSI2/719 


3600 


SuSI2/717 


1800 


NOT observations 


















IRAS 08339+6517 


ALFOSC/70 


2400 


ALFOSC/78 


3000 


ALFOSC/113 


5000 


ALFOSC/17 


3300 


NGC 6090 


ALFOSC/53 


4200 


ALFOSC/78 


2400 


ALFOSC/40 


3600 


ALFOSC/17 


2100 


Table 2. Ground-based observations of our six targets. Northern targets have been observed with the Nordic Optical Telescope 



and southern ones with ESO New Technology Telescope. Ha, H/3 online and respective continuum observations are listed with 
the instrument and filter name and the exposure time (in seconds) in each band. 



Section |2] we describe the observations and the data reduction, 
in Section[3]we present the results, Section|4]is dedicated to the 
discussion of these results and finally in Section [5] we present 
our conclusions. We assume throughout this paper a cosmology 
of //o=72 km s-i Mpc"', Qm=0.3 and ^a^O.V. 



2. OBSERVATIONS AND DATA REDUCTION 

Our sample contains six local starburst galaxies that have been 
hand-picked in order to cover a range of relevant intrinsic param- 
eters. Based upon UV-properties, the selection covers a range 
of dust content, luminosity and the variety of observed Lya 
profiles: firstly, the sample consists of four Lya emitters from 
|Calzetti & Ki nney (1992), among which candidates represent- 
ing a range in line profiles are selected ( Kunth et al.|1998 [Mas- 
Hesse et al.|2003 1. Two known Lya absorbers are also included 
from |Thuan & Izotov (|1997| l. General information and proper- 



ties are given in Table 



2.1. HST observations 

Two general observer programmes, with the Hubble Space 
Telescope (HST), have been devoted for the observations of the 
six targets: GO9470 that makes use of the Solar Blind Channel 
(SBC) of the Advanced Camera for Surveys (ACS) in order 
to observe Lya online (F122M) and continuum (F140LP), and 
GO 10575 using the High Resolution Camera (HRC) and Wide 
Field Camera (WFC) to obtain Ha, near-UV and optical con- 



tinuum broad band observ ations. Detailed des cription of these 
observations is presented in Osthn et al. ( 2008 1. 

All images are 'drizzled' using the MULTIDRIZZLE task 
in STSDAS package under NOAO/IRAF onto the same pixel 
scale (0.025" pix"') and the same orientation. Remaining shifts 
are corrected with GEOMAP and GEOTRAN, and cosmic rays 
removed using CREDIT task. The images for each target are 
then convolved to the same Point Spread Function (PSF) using 
DIGIPHOT/DAOPHOT, with the worse PSF image as reference. 



The production of continuum subtracted Lya images from 
F122M (online) and F140LP (offline) is not trivial and needs 
sophisticated techniques for many reasons. The effective wave- 
length of the continuum filter is rather far from the online filter 
(AA/A ^ 0.22) and the UV continuum between the two filters de- 
viates significantly from a power law (/) oc A^) and is sensitive to 
age and E(B-V). Neglecting these considerations and describing 
the behavior of the UV continuum near Lya with a power law 
yielded results ranging from absorption to strong emission ac- 
cording to the values set to /3. Preliminary results, bylKunth et al. 
(|2003), pointed out these limitations and the need for an elabo- 
rated subtraction method. Subsequently, [Hayes et al.| ( [2005| l and, 
in more details, Hayes et al. ( 20081) presented a reliable extrap- 



olation from F140LP to F120M described by the Continuum 
Throughput Normalization factor (CTN). To estimate the CTN 
in each pixel, all the images from F122M to F814W are used to 
fit the Starburst99 spectral evolutionary models (Leitherer et al. 
1999 Vazquez & Leitherer|2005 1. Summarily, the method uses 
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filters that sample the UV continuum slope and the 4000A break 
to fit burst age and stellar E(B-V) using minimisation. Then 
the nebular contribution and the stellar components are treated 
independently. Indeed, Ha image is used to constrain the nebu- 
lar gas contribution to the total SED, while V and I images allow 
to estimate the stellar part. The gas spectrum is subtracted and 
age and mass are fitted in two stellar components. 

To generate Starburst99 models, we set the metallicity as de- 
rived from observations of each object: for IRAS 08339+6517 
and NGC 6090 Z = 0.02, for Haro 11 and ESQ 338-04 
Z = 0.004, and for SBS 335-052 and Tololo 65 Z = 0.001. 
Though Hayes et al. (2008) show that an error below 50% on 
the metallicity does not significantly affect the continuum sub- 
traction, the metallicity estimate remains the factor driving the 
accuracy of the subtraction method. 

A standard [Saipeter| ^T955\ IMF is used (a = 2.35, dN = 
M'"dM) in the range O.IMq to 120Mo. Multi-stellar component 
fitting, used for continuum subtraction, is almost completely in- 
sensitive to the IMF slope and mass range. Finally, an instan- 
taneous burst scenario is adopted as it is more appropriate for 
individual pixels but the choice of a constant star formation sce- 
nario would not affect global photometry of the galaxies or the 
contribution of the underlying stellar population. 

2.2. ground-based observations 

The southern targets of our sample (cf. Table [2]) have been 
observed using the New Technology Telescope (NTT) at La 
Silla (ESQ) during the nights of 18, 19 and 20 September 2004 
(apart from Tololo 65 which was observed in service mode on 
28 January 2003). Narrow-band imaging was performed for all 
targets in Ha, H/3 and [Oiii]/15007A, and nearby continuum for 
each line. The first night of the run, 18th Sept, the seeing was 
consistently below 1.2", yet the presence of thin cirrus prevents 
a direct calibration using standard stars. Observational condi- 
tions on the night of 19 Sept were excellent: photometric sky 
and sub-arcsec seeing. On the last night the seeing was worse, 
exceeding 2", although the photometric quality was still good. 
Spectrophotometric standard stars FeigellO, LDS749B, GD50 
and GD108, selected from ,Oke, (1990j catalog, were observed 
at regular intervals during each night and in each filter. Both the 
ESO Multi-Mode Instrument (E MMI) (|Dekker et al.|l 986) 1 and 
Super Seeing Imager 2 (SuSI2) P'Odorico et al.J|199^) were 
used interchangeably. A binning of 2x2 pixels was used, giving 
plate-scales of 0.332" pix"' and 0.161" pix"', and fields-of-view 
of 9.1 X 9.9' and 5.5 x 5.5' for EMMI-R and SuSI2 respectively. 
NTT observations are summarized in table [2] together with ESO 
filter codes. The good seeing observations of the night 18 Sept 
are calibrated using secondary standard stars in the field from 
the photometric night 19 and/or 20 Sept. 

The remaining northern targets (IRAS 08339+6517 and 
NGC 6090) were observed with the Nordic Optical Telescope 
(NOT) at La Palma during the nights of 14, 15 and 16 
February 2004, using the Andalucia Faint Object Spectrograph 
and Camera (ALFOSC). Unfortunately, the observational condi- 
tions were not photometric during the observing run. For NGC 
6090, we used observations taken on March 28th 2006 under 
photometric conditions. The calibration was performed using 
BD33 standard star. Our science images were then calibrated us- 
ing field stars from 2006 observations in each filter. For IRAS 
08339+6517 , in the same purpose, we made use of HST im- 
ages. Firstly, calibrated Ha images are rescaled to the ALFOSC 
pixel-plate (0.19" pix '). Aperture photometry is then obtained 



and used for the calibration of the ground-based Ho- observa- 
tions. The H/3 continuum images are calibrated using rescaled 
F435 W images and Hy6 online images are then calibrated with re- 
spect to the continuum images using photometry of several stars 
in the field of view. Global estimates for IRAS 08339+6517 ai-e 
the least accurate of the galaxy sample due to the less-than-ideal 
calibration of H/3. However, pixel-to-pixel measurements suffer 
from similar uncertainties, therefore pixel-to-pixel trends are as 
reliable as they are in any other target, irrespective whether the 
exact E(B-V) is perfect. 

Data were first reduced using the standard NOAO/IRAF proce- 
dures. Images were bias-subtracted, and flat-field corrected us- 
ing well exposed sky and dome flats. Images were registered 
using field stars and the GEOMAP and GEOTRAN tasks, and 
images from the same bands were coadded with inverse vari- 
ance weighting. Finally, all reduced frames for each target were 
smoothed to the PSF of the worst seeing image using GAUSS. 

Hff images were corrected for [Nn] contamination and both 
Ha and H/3 images corrected for underlying stellar absorption. 
[Nil] contamination is estimated from Nn A65^3A/Ha found in 
Hterature (see Table|3]l and the ratio Nii /i6584A/Nii /16548A ^ 3 
( |Osterbrockri9 89). 

Stellar absorption is estimated using the best-fitting 
Starburst99 spectrum given by the procedure described in Sec. 
|2.1| Latest stellar atmospheres models are used in Starburst99 
stellar libraries (Martins et al., 2005), including full line- 
blanketing, which allows the best estimate of equivalent width 
of the Ha and H/3 stellar absorption features. 

Finally, using the Balmer decrement Ha/H/3, assuming a case 
B intrinsic line ratio of 2.86 ( |Osterbrock| 1 989[ ) at a temperature 
of 10,000 K, we can derive the E(B-V) in the gas phase using 
the parameterization of Cardelli et al. ( 1989| l and the equation: 



E(B - V)Ha/HI3 



2.5 X log(2.86/Rob,) 



(1) 



where Robs - fual fH/} is the absolute observed flux ratio, and 
k(/lQ.), k(/la) are the extinction curves at Ha and Hf3 wavelengths 
respectively. According to |Cardelh et al.| ( [1989] ), k(/l„) ~ 2.63 
and k(^^) -3.71. 



Target 


[Nil] ,I6583A 
Htr 


Reference 








Haro 1 1 


0.189 


'Bergvall & Ostlin'('2002k 


ESO 338-04 


0.031 


Bergvall&Ostlin (20021 


NGC 6090 


0.411 


Moustakas & Kennicutf]'(|2006 


IRAS 08339+6517 


0.25 


Margon et al. ( 1988j) 


SBS 0335-052 


0.003 


izotov et al. 1 1997| 


Tololo 65 


0.005 


izotov et al. 1^20041 



Table 3. Ha correction for Nii contamination with references 
where the ratio has been taken from. 



2.3. Uncertainties 

The calibration of individual images in Ha and Hf) lines are 
subject to classical errors in the data reduction procedure. The 
origin of the uncertainty comes from the flux measurement of 
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the calibration standard stars. This accounts for about 5% in our 
observations. Errors of the order of few percent may be produced 
by other effects, such as residuals in the flat-field corrections, 
residuals in the sky background and continuum subtraction. We 
make the assumption that all these effects will lead to an eiTor on 
the flux measurement of the order of 10%. Moreover, we are pro- 
ducing extinction maps by dividing Ha by li/3 images, and, for 
some galaxies, these images are obtained with different instru- 
ments (i.e. different plate-scales and orientations), the errors on 
images alignment and registration may lead to quite significant 
uncertainties, since we aim, in our science analysis, to investi- 
gate the Lya and dust amount variations on a pixel scale. The 
procedure to estimate the impact of such misalignments was the 
following: 

The RMS of the alignment fit given by GEOMAP is about 
0.2 pixel. Consequently, to estimate the misalignment (and only 
to that purpose), we have rebinned Ha and H/3 images to a new 
pixel size (1/5 of the original pixel size). Then we have created 
new images with artificial shifts by -h1,0 or -1 pixel in x and 
y directions. This gives us 9-image data cube for both Ha and 
H/3. Using all combinations of Ha and Hfi shifted images we 
derived the ratio Ha/Hfi and constructed an E(B-V) data cube. 
We eventually computed the standard deviation of the extinction, 
C£(s-v), at each pixel [x,y]. 

Misalignment errors are also a concern for the Lya image 
with respect to Ha and online images with respect to contin- 
uum ones. This could result in errors on the Lya/Ha ratio and 
the equivalent widths. The estimation of these quantities takes 
into account this uncertainty. Besides, we have accounted for 
the uncertainties of the Lyo' continuum subtraction procedure by 
assuming a statistical error of 10% of the Lyo' flux. This uncer- 
tainty depends, in reality, on the signal-to-noise in each resolu- 
tion element where the continuum is fitted, and the value adopted 
here corresponds to S/N=5%. The final eiTor bars are shown on 
the scatter plots in Sect. [3] For the sake of readability we plot the 
mean eiTors instead of individual eiTors for each point. When the 
emission is dominated by the conglomerate points of the diffuse 
halo (which is the case for most of the plots), we show rather 
the mean computed between only the maximum and minimum 
values, in order to cover the full range of eiTor variations. 



3. Analysis 

HARO 11 — The Ha image in Fig. [Tjin the first column shows 
a complex morphology with three main star-forming condensa- 
tions (Kunth et al. 2003). The continuum subtracted Lya image 
does not delineate this morphology, showing Lya in emission in 
only knot C, whereas it is seen in absorption in knot A and B. 
The decoupling of Lya from continuum is clearly observable in 
the bottom frame, which represents Lya map at HST resolution 
overlayed by FUV (1500A) contours. The emission exhibits two 
different components consisting of a central bright knot and a 
low surface brightness diffuse emission. It seems that, examin- 
ing the extinction map, the diffuse component is not quite regu- 
lated by the dust amount. Moreover, the bright Lya emission in 
knot C, corresponds to a high extinction region. 

This galaxy is a well known Lya emitter ( Kunth et al.|1998 , 
while the detection of Lyman continuum leakage by Bergvall 
let al.| ( |2006| l is still controver sial ([Grimes et al.J|2007| l. It has 
been studied in more details in [Hayes et al.|(|200T 



ESO 338-04 — In the second column of Fig. [T] the Lya image 
shows three main absorption regions and a surrounding bright 



emission. The absorption sites coiTespond to relatively dusty re- 
gions of the galaxy seen in the extinction map, which traps Lya 
photons, while the emission is not correlated with the dust con- 
tent. The last component is, again, the diffuse emission around 
and overall the galaxy, with a low surface brightness, which is 
showing up the resonant decoupling of Lya photons. Many dust 
features are clearly visible in the E(B-V) map, with a clumpy- 
like structure, and following roughly the Ha structure. 

The Lya image is produced by matching the HST/ ACS im- 
age to the NTT resolution. This results in spreading out the cen- 
tral absorption to more extended region, revealing small new ab- 
sorption features around and dimming the surface brightness of 
the emission component. The bottom frame shows again how 
Lya is uncorrected with the FUV continuum which traces un- 
obscured star formation sites. 



NGC 6090 — Third column in Fig. [T]shows that the interacting 
system NGC 6090 exhibits Lya emission from each component, 
at a distance of about 6" from each other. The emission is peak- 
ing around low extinction regions, and the overlapping region 
between the two components, where no emission is seen appears 
very dusty. The extinction map of NGC 6090 shows dust pattern 
like a spiral structure and no evident correlation with the ionised 
gas traced by the Ha emission is seen. 

The main Ha structures correspond to the Lya emission 
components, although the largest Lya emission represents a 
small region in Ha, and vice-versa. This discrepancy may be 
due the large amount of dust in the upper left component that 
could destroy a part of Lya photons. Knot A appears also dustier 
{EB-v,gds ~ 0.75) than knot B (£B-v,gas ~ 0.55). 

IRAS 08339+6517 - This nuclear starburst shows a spiral 
structure conspicuous on the Ha (Fig [2]i and FUV continuum 
images (the latter is not shown here). However the Lya image 
does not resemble any of them, showing a central bright emis- 
sion and the ubiquitous halo component. The FUV contours ex- 
hibit much more details in the arms of the galaxy consisting in 
many star clusters where Lya is absent. 

The extinction map indicates a dust-free central spot. The 
dust distribution has no clear relationship with the emission 
maps. 

SBS 335-052 — The dust distribution is apparently correlated 
with the Ha emission in this galaxy. In Fig. |2] the brightest re- 
gion in Ha presents the most important amount of dust. Are also 
visible in both images a shell structure and a dust-free region 
toward the S-E of the bright region. 

The bright dusty spot coiTesponds to a relatively high Lya 
absorption region. The galaxy shows Lya only in absorption sur- 
rounded by a faint diffuse emission. It appears that, despite a 
bright Ha or FUV continuum emission, no Lya photons manage 
to escape directly without scattering on neutral hydrogen. 

Tololo 65 — In the last column of Fig. |2] no Lya structure 
can be seen for this galaxy, though the HST image shows some 
emission features (bottom frame), they have been smoothed to 
the NTT resolution. Degrading the resolution of the HST im- 
age leads to ghost features since the emission becomes as weak 
as the background level and is swallowed up. These artefacts in 
the background, although less dramatical, are also present in the 
other galaxies. The Ha map presents a two-components emis- 
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Fig. 1. Galaxy imaging (Part 1): From top to bottom: Ha Lya, E(B-V)jfl5 map, and Lya as seen by HST overlayed with FUV 
contours. Inverted logarithmic scale is used, showing emission in black and absorption in white. The extinction map is overlayed 
with the mask generated following the description in Sect. |4] Regions out of the contour are excluded from our study. The dustiest 
regions are in black. /rom left to right with the FoV and spatial scale in parentheses: Ground-based images (top three rows): Haro 
11 (17", 0.4 kpc/"), ESO 338-04 (17", 0.2 kpc/"), and NGC 6090 (19", 0.57 kpc/"); HST images (last row): Haro 11 (13"), ESO 
338-04 (15"), and NGC 6090 (21"). North is down and East to the right. 




Fig. 2. Galaxy imaging (Part 2): From top to bottom: Ha, Lya, E(B-V)(,as map, and Lya as seen by HST overlayed with FUV 
contours. Inverted logarithmic scale is used, showing emission in black and absorption in white. The extinction map is overlayed 
with the mask generated following the description in Sect. |4] Regions out of the contour are excluded from our study. The dustiest 
regions are in black. /rom left to right with the FoV in parentheses: Ground-based images (top three rows): IRAS 08339+6517 (15", 
0.38 kpc/"), SBS 335-052 (4", 0.27 kpc/"), Tololo 65 (8", 0.18 kpc/"); HST images (last row): IRAS 08339+6517 (14"), SBS 
335-052 (5"), Tololo 65 (8"). North is up and East to the left. 
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Fig. 3. Haro 1 1 scatter plots. Top-left: Lya surface brightness as a function of the extinction determined from the Balmer decrement. 
Regions of knots A, B, C and the diffuse emission are marked in the figure and represented with different colors and different 
symbols: Emission is in red, absorption is in blue (knot A) and green (knot B), and the diffuse emission component is in cyan. For 
the remaining plots, the same color code is used for the respective regions. Top-right: Lya equivalent width E(B-V). Bottom-left: 
Ha equivalent width vs E(B-V). Bottom-right: LyafHa ratio vs E(B-V). Overplotted on all figures the error bars corresponding to 
the uncertainties estimated in Sec.P 



sion that resembles the structure of the dust content in the E(B- 
V) map. 

In order to investigate the Lye emission variations and pos- 
sible correlations with dust or other parameters at the smallest 
possible scale, we have produced scatter plots from the images. 
In extinction maps obtained in Fig.[T]each point represents one 
pixel in the galaxy region which has been isolated by masking 
the background at 5cr level. Regions of interest, such as star- 
forming knots and emission or absorption features, are picked 
up using circular apertures and represented with different colors 
and symbols on the figures. 

3. 1 . Blended emission and absorption systems 

Haro 11 - The first plot in Fig. [3]presents the correlation be- 
tween the Lya emission and the extinction determined from the 
Balmer decrement tracing the dust in the gas phase (EB-v.gRs)- 
The color-code represents different regions of interest consist- 
ing of circular apertures centered on the three main knots of the 
galaxy, which are marked on the Lyo' image. We can see a dif- 
fuse emission component extending up to Eb-v ~ 1-5. Knot C 
shows a bright and spread emission with a mean extinction of 



0.48, while the absorption is essentially localised around knot A 
and B. The mean extinctions in each knot are derived using the 
ratio of integrated Ha and HyS fluxes within the corresponding 
apertures. 

The presence of two distinct emission components is indica- 
tive of two different physical processes implied in the escape 
of Lya photons. Firstly the diffuse component that shows the 
resonant decoupling of Lya photons resonantly scattered on hy- 
drogen atoms until they escape far away from there production 
sites (thus, experiencing a large range of extinction). On the 
other hand, the emission from knot C is less spread and rep- 
resents photons escaping directly from this small region with a 
mean extinction Eb-v ~ 0.48. The Lya resonant decoupling is 
also visible in Fig. |4] representing the correlation between the 
Lya and the Ha fluxes in log-scale, hence, showing only pos- 
itive pixels (in emission). We observe a first component at a 
low and roughly constant Lya surface brightness (around 10"'"^ 
erg s ' cm"^ arcsec"^), independent of Ha emission. Thanks to 
their resonant scattering, Lya photons reach regions where non- 
resonant photons, such as Ha, are absent, making the Lya/Ha 
ratio higher than the Case B level (represented by a dashed line 
in the figure). The second component at higher Lya and Ha 
fluxes (/Lyo- > 10"''* erg s"' cm"^ arcsec"^ ) is always below the 
predicted recombination value. These pixels represent regions 
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Fig. 5. ESO 338-04 scatter plots. Top-left: Lya surface brightness vs E(B-V). Emission from the central region is in red. Absorption 
in the central region (knot A) is represented in blue (knot A), the surrounding emission is in red, and the diffuse emission component 
is in cyan. Same color code applied in all the plots. Top-right: Lya equivalent width vs E(B-V). Bottom-left: Ha equivalent width vs 
E(B-V). Bottom-right: LyafHa ratio vs E(B-V). Error bars correspond to the uncertainties estimated in Sect. 2.3 
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Fig. 4. Pixel distribution of the \jya against Hot surface bright- 
ness. A logarithmic scale is used and hence shows only positive 
contribution (emission). The dashed line represents the case B 
recombination ratio. The galaxy name is given in each plot. 

where Lyor photons escape directly from their production site 
where Ha is also produced. 



That we see Lya in emission from knot C with Zsg-v ~ 0.48 
whereas absorption is seen in knots A and B with E^^y ~ 0.2 and 
0.41 respectively is interesting. Indeed Lya photons manage to 
escape from regions with higher extinction than those of a pure 
absorption. The dust content is clearly not the main driver in the 
escape process of Lya photons. A great covering of static Hi 
column density in knot A and an expanding neutral ISM and/or 
ionised Hii holes in knot C may lead to this observation. 

Looking at Lya equivalent width, we observe that the diffuse 
emission shows relatively high EWi,ya whereas the pure emis- 
sion equivalent width in knot C (in red) is much weaker. This 
suggests that a hard FUV radiation could create ionised holes 
through which Lya photons may escape, in a inhomogeneous 
distribution of Hi and dust. In this case of multi-phase ISM, it 
has been shown (Neufeld,1991, ,Hansen & Peng Oh_2006 ) that, 
thanks to their scattering on the dusty Hi clumps, Lya escapes in 
an easier way than non-resonant photons. We also observe that 
the diffuse emission (in cyan) corresponds to the highest equiv- 
alent width observed (EWiya higher than 200 A), since it repre- 
sents photons that have scattered far away from their production 
sites and escape where the Lya continuum is lower. This decline 
in the emission is, again, symptomatic of the resonant nature of 
Lya photons. Indeed, when we plot the equivalent width of Ha 
against the extinction, we do not see any correlation, as we ex- 
pect for non-resonant lines, since the online and the continuum 
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Fig. 6. IRAS 08339+6517 scatter plots. Top-left: hya surface brightness vs E(B-V). The nuclear Lya emission component is rep- 
resented in red, and the diffuse emission component is in cyan. Same color code applied in all the plots. Top-right: Lya equivalent 
width vs E(B-V). Bottom-left: Ha equivalent width vs E(B-V). Bottom-right: Lya/Ha ratio vs E(B-V). Error bars correspond to the 
uncertainties estimated in Sec. 12.31 



photons are regulated by the same physical processes according 
to the dust content. The last figure shows how the Lya/Ha ratio 
evolves according to the amount of dust. In a classical view, con- 
sidering only the selective extinction at the two wavelengths, and 
a case B intrinsic ratio of 8.7 ( B rocklehurst| 1 97 l] l, we expect to 
have an exponential decline represented by the dark curve. The 
resonant nature of the Lya photons leads to a different result. We 
observe a high dispersion for the halo component and an emis- 
sion from knot C above the predicted level at higher extinction, 
which might support the view of a scattering in inhomogeneous 
ionised ISM that favors preferentially the escape of Lya photons. 



ESO 338-04 - This galaxy shows a similar pixel distribution 
in Ly o-Eb-v space to that of Haro 1 1 (Fig.|5]l. We see both emis- 
sion and pure absorption features at, approximately, the same ex- 
tinction (mean Eg-v ~ 0.22 and 0.23 for absorption and global 
emission respectively). The different regions are marked in the 
scatter plots with different colors consisting of a central absorp- 
tion (in blue) surrounded by emission features (red). We also ob- 
serve a halo of diffuse emission surrounding the starburst regions 
and independent of the extinction. However the diffuse emission 
does not exceed 5 x lO"'** erg s"' cm"^ arcsec"^, it accounts 
for about 70% of the total Lya emission from the galaxy. The 
same value was found by Hayes et al. ( 2005 1, although a differ- 
ent masking was used. 



The absorption is seen only in the knot A (following the 
nomenclature of Hayes et al. 2005 1 superimposed on a region 
with an average extinction of Eb-v ~ 0.2 (mean extinction cal- 
culated, as for Haro 1 L using integrated Ha and Hj3 fluxes over 
this knot). 



Ostlin et al. 



(2003) from longslit spectroscopy mea- 
surement of Ha/HyS (the slit being positioned on the East- West 
direction) found E(B-V) values within a range comparable to 
our result (Eb-vshs ~ - 0.25), except toward the east edge of 
the galaxy where a dustier region is seen in our extinction map 
{EB-v,gas ~ 0.4). Same extinction as in knot A is observed in its 
surrounding regions, where we found Lya in emission. There 
is a trend of decline in the emission (though, with high disper- 
sion) (until Eb-v ~ 0.6). The radiative transfer of Lya photons 
in a static (or nearly static) ISM may lead to the present situa- 
tion: In the central region (knot A), Lya photons are produced 
and immediately absorbed by neutral hydrogen. They are re- 
emitted according to space and frequency redistribution prob- 
ability. Eventually, Lya photons resonantly scatter in the wings 
until they reach a frequency which is far enough from the line 
center (where the absorption probability is close to unity) to be 
able to escape from the neutral medium. This also leads to a 
diffusion in space which explains the emission seen around the 
central absorption at the same extinction. Actually, this effect 
is also seen in Haro 11 with the HST resolution (Fig. [T]i, but 
does not appear at the NTT resolution due to the smoothing ef- 
fects. It seems that, as shown by Lya radiative transfer models 
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Fig. 7. NGC 6090 scatter plots. Top-left: Lya surface brightness vs E(B-V). Emission from knot A is in red, and from knot B is in 
green. Diffuse emission component is represented as usual in cyan. Same color code applied in all the plots. Top-right: Lya equiv- 
alent width vs E(B-V). Bottom-left: Ha equivalent width vs E(B-V). Bottom-right: LyafHa ratio vs E(B-V). Error bars correspond 
to the uncertainties estimated in Sec. 12.31 



(e.g. Verhamme et al.|[2006) l, Lya photons do not easily escape 
from where they have been produced but scatter in space and 
frequency avoiding absorption by HI atoms at line center (vq - 
2.46 X lO'-'^ Hz). 

We note in Fig. |5] very high Lyor equivalent width, ranging 
from -10 A for the absorption (blue) to ~250 A in the surround- 
ing emission (red). Such high EW might be due either to the 
resonant scattering mechanism that allows Lyo' photons to travel 
where continuum photons can not, or to a multi-phase configu- 
ration of the ISM where, [Hansen & Peng Oh ( 2006[ ) noted for 
reasonable HI column density and dust amount, that continuum 
is preferentially extinguished, boosting the initial EW easily by 
a factor of 2-5. This configuration may also explain the high 
EWi^yo observed in some high-redshift galaxies, extending up 



to 150 A in LALA z ~ 5.7 sources ( |Rhoads et aL]|2003[ ) in the 



case of spectroscopically confirmed candidates (higher possible 
equivalent widths were found in the imaging survey). On the 
other hand. Ha equivalent width map shows similar distribution 
to that of Haro 1 1 without any correlation. 

The last frame shows the evolution of Lya/Ha ratio with the 
extinction. In the central region around knot A, Lya/Ha follows 
loosely the theoretical curve (marked in black). 

However, the absorption region does not show any trend of 
decrease in LyafHa ratio. Rather, we see a dispersion, as seen in 
the Lya-Es-v plot, ranging from Eb-v ~ to 0.5. 



3.2. Emission systems: NGC 6090 & IRAS 08339-h6517 

Both of these galaxies exhibit Lyo- only in emission, with little 
signs of direct absorption, as we have seen on Fig. [T] and [2] As 
usual, the first plot in Fig. 6] and |7] shows Lya flux function of 
the color excess Eg-v- In both systems, we do see the diffuse 
component independent of Eg-v at low Lya flux. 

Although, the calibration procedure used for IRAS 



08339-1-6517 data is not accurate (see section 2.2 1, the average 
extinction found for this galaxy ( Eg-v ~ 0.12 ) is still close to 
the value estimated by [Gonzalez Delgado et al. 



T998] l. 



STIS spectroscopy (Ma s-Hesse et al.[[2003 i revealed a P- 
Cygni profile, with a red wing shifted by ~ 300 km s ' with 
respect to the Hii region velocity determined from the optical 
emission lines. The extension of the neutral gas shell is found to 
be very large, with a diameter of around 10 kpc, which is larger 
than the aperture used here to isolate the emission region in the 
center of the galaxy (~ 2 kpc). Thus, the central emission (red 
component) is observed through an expanding shell approaching 
us at a velocity of ~ 300 km s"' . 

Same Lya profile has been observed for NGC 6090 by 
Gonzalez Delgado et al. ( 1998[ l, with a velocity offset with re- 
spect to blueshifted interstellar absorption lines, indicative of 
large scale high-velocity outflows of gas, of several hundred 
km s" ' . We differentiate two emission blobs (red for A and green 
for B in Fig.[7]i. Component A appears at a mean E(B-V) around 
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Fig. 8. SBS 335-052 scatter plots. Top-left: Lya surface brightness vs E(B-V). The damped central absorption is represented in red, 
and the diffuse emission component is in cyan. Same color code applied in all the plots. Top-right: Lya equivalent width vs E(B- 
V). Bottom-left: Ha equivalent width vs E(B-V). Bottom- right: Lya/Ha ratio vs E(B-V). Error bars correspond to the uncertainties 
described in Sec. 12.31 



0.74 whereas the second one atE(B-V) ~ 0.54. The diffuse emis- 
sion, also in this galaxy, is the only component that attains the 
recombination level (Figure |4]i while its bulk stays quite below. 

Again, regarding these two galaxies, the ISM kinematics 
must play an important role in the escape of Lya photons through 
neutral gas shells, making Lya less sensitive to the dust con- 
tent. It leads to the observed dispersion in the emission over a 
large range of extinction for these two galaxies. To confirm this 
visual inspection, we have performed a Spearman's statistical 
test which gives the probability that a correlation between Lya 
and Eb-v exists. The Spearman's correlation coefficient p can 
take values from -1 to -i-l. A value of +1 shows that the vari- 
ables are perfectly correlated with a positive slope, a value of -1 
shows that the variables are anti-correlated (negative slope), and 
a value of shows that the variables are completely indepen- 
dent. While we would expect an anti-correlation between these 
two variables, we found p ~ 0. 15 for both galaxies (whereas p ~ 
-0.2 and -0.05 for Haro 1 1 and ESO 338-04 respectively), con- 
firming that Lya is less sensitive to dust in these systems. For 
SBS 335-052 we found p ~ -0.56. 

The Lya equivalent width appears to be not correlated with 
the dust in IRAS 08339-H6517 . This absence of correlation is 
also seen when observing the evolution of the line ratio Lya/Ha 
which is around 2 and stays constant as the extinction increases 
supporting the importance of the ISM kinematics. 



The situation is slightly different for NGC 6090, where in 
knot A we see values more clearly above the theoretical Lya/Ha 
curve (as seen for Haro 11). In addition to kinematics, we may 
be in presence of clumpy ISM configuration as well (as it is sug- 
gested by the images), that allows Lya to get out through the 
ionised inter-clouds medium. Furthermore, the two knots show 
different equivalent widths: EWiya ~ 54 A in knot A and 20 A 
in knot B. Note that, the diffuse emission presents higher values 
because of the numerous resonant scatterings that have experi- 
enced these photons. This is still valid for Lya/Ha ratio. 



3.3. Damped absorption system: SBS 335-052 



HST/GHRS spectroscopy ( Thuan et al.|1997| l has revealed these 
metal-deficient BCDs to be damped Lya absorbers. They don't 
show, indeed in our images, any direct emission but a low sur- 
face brightness diffuse halo. The Lya photons manage to escape 
from these galaxies after multiple scattering events. Figure |8] 
shows how Lya is related to different physical parameters. We 
see, for SBS 335-052, in the first frame, a weak diffuse compo- 
nent which accounts for the whole Lya emission. In a classical 
vision, where dust is the main parameter responsible of the de- 
struction of Lya photons, we expect a declining relationship be- 
tween Lya flux and the dust extinction. It is precisely what we 
observe in this galaxy. A strong absorption (fiya 3 x 10"^^ 
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erg s ' cm"^ arcsec"^) is seen, with a weak declining trend in the 
range i Eb-v i 0.7. This result suggests that the dust is play- 
ing, in this case, an important role in the escape of Lya photons. 
Very Large Array 21 cm observations (Thua n & Izotov|[l997| l 
show the BCD to be embedded in a large Hi cloud. The Hi col- 
umn density in the GHRS aperture derived by thes e authors is 
very large, A^(Hi) = 7.0 x 10^' cm"-. According to 



Mas-Hesse 



et al. ( 2003| l and the evolutionary models of Tenorio-Tagle et al. 



( 1999 1, this galaxy is a very young starburst of which the age 
of the stellar population is too small to have ionised the whole 
surrounding medium, leading to a great amount of neutral gas 
covering the massive stars. In addition, the Hi cloud is static with 
respect to the central Hii region as the 21 cm and the emission 
line velocities are in good agreement. As a result, the Lya pho- 
tons after a large scattering process on hydrogen atoms increase 
their mean path and the probability to be absorbed by dust grains. 
The combination of a great Hi column density in front of the 
Lya production sites and the absence of kinematics in this neu- 
tral envelope makes Lya photons very likely to be destroyed by 
dust, hence turning the dust extinction into a very important Lya 
escape regulator. This configuration, with the related dust corre- 
lations, was observed for the knot A in Haro 11 . It also appears 
that the Lya equivalent width is declining with the extinction. 

TOLOLO 65, in the same way, shows very weak Lya 
emission through a diffuse component without direct emission. 
GHRS spectrum (Thuan & Izotov 1997) shows, with an accept- 
able signal-to-noise ratio although poorer than that of SBS 335- 
052, that this galaxy is a pure Lya absorber. This is consistent 
with the very faint and diffuse emission found here, which is still 
not significantly greater than the background level where the de- 
graded resolution w ipes out the absorp tion seen in the original 
HST/ACS image by|6sthn et al](|2008|. 



4. Discussion 

Rather than the detailed description performed in the last section, 
for each galaxy to investigate how Lya emission is related to 
other parameters, especially dust, at a small scale; we discuss 
hereinafter the global characteristics of the galaxies based upon 
the integrated quantities in defined apertures. These apertures 
are defined for all the galaxies by masking regions that shows 
f(^8) below a threshold of 5 times the standard deviation of the 
background in the Hj6 image. This study of course does not allow 
any meaningful statistics with only six galaxies in our hands, 
that have been however selected in order to sample the largest 
possible starburst physical parameters. Yet we can outline some 
interesting trends when compared to the small-scale approach 
and how far some properties could be smoothed, or not, at the 
galaxy scale. 

4.1. The role of the dust in Lya obscuration 

We show in Table |4] the photometric properties of our sample, 
integrated in an aperture defined by a mask based upon a H/3 flux 
threshold above the sky level, in order to get the extended emis- 
sion of the Balmer lines and the Lya diffuse emission. As men- 
tioned before, all the quantities presented have been corrected 
for galactic extinction using Schlegel et al.| ( |T998[ l method. The 
integration of the Lya flux within the mask aperture results in 
five candidates of our sample being emitters and only one net 
absorber (SBS 335-052). Defining this galaxy as a net absorber 
means that the sum of the flux overall the galaxy gives a negative 
result, however we do see emission in some regions. Though, the 
flux and the measurements derived from are very sensitive to the 
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Fig. 9. Scatter plot of integrated quantities, representing in each 
galaxy the Lya equivalent width as a function of the extinction. 
E(B-V) is derived from the ratio between the integrated Ha and 
HyS fluxes. The legend of symbols assigned to the galaxies are 
showed in the inset. We have included in this plot the two com- 
ponents of NGC 6090 that we have considered separately. 



aperture size. We can expect that we are losing a part of the weak 
diffuse Lya emission and the most extended Ha or Hfi emission, 
which could be attained by deeper observations. Thus, measure- 
ments such as the escape fraction of Lya may change according 
to the adopted mask size since Lya can scatter further from the 
production sites than Balmer or continuum photons. 

The Lya/Ha ratio ranges from -1.12 to 1.36, showing Lya 
emission much weaker for all the sample than predicted by re- 
combination theory even, in most cases, when corrected for the 
differential extinction at these different wavelengths. Previous 
observations (Ter levich et al.||1993t [Giavalisco et al.||1996| l 
yielded the same conclusions. The dust is just the final stage 
of the process responsible for the obscuration of Lya photons, 
after the resonant scattering in an homogeneous medium that 
increase their mean path, which makes the recombination ratio 
not only regulated by the dust (as seen from the Ha/HyS ratio 
which traces the nebular dust). Setting the Lya equivalent width 
against the different parameters of Table |4]ends up with no clear 
trends. Plotting the Lya equivalent width against E(B-V) allows 
us to probe the difference in extinction between resonant and 
non-resonant radiations, knowing that EWiya is not affected by 
the selective extinction (i.e. independent of the dust extinction 
curve). We see, indeed, in Fig. |9j a rather scattered data points 
and no well-observed correlation. Since we are dealing with res- 
onant radiation investigated from only one line-of-sight one ex- 
pect that the geometry and the distribution of the dust layers 
around the emitting regions may affect the observed scatter of 
EWiya according to the extinction. This relation depends also 
on the intrinsic EWs. 

When we calculate the dereddened ratio (Lya/Ha)c, rather 
than EWijya, in Table HI we observe that it lies below the re- 
combination value (8. 777 expected if the extinction is only due to 
dust, for five galaxies. Only NGC 6090 shows (Lya/Ha)c that 
exceeds this recombination level. Since the resonant Lya radi- 
ation is spatially decoupled from continuum or Balmer lines, 
the Lya photons may suffer different extinction than traced by 
the Balmer decrement and lead to this overestimation of the 
reddening correction. In other words, the Lya photons that are 
detected and survived to numerous scattering and various ex- 
tinctions, have been travelling through ISM regions where dust 
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Target Aperture fLy„ fua fiys Lya/Ha (Lyci'/HQ')f lla/H/3 EW(Lya) EW(Ha) 

Size (erg s"' cm"~) (erg s"' cm"^) (erg s"' cm"^) (A) (A) 



Haro 11 


143 (23) 


1.3e-12 


2.3e-12 


5.48e-13 


0.57 


5.48 


4.14 


22.8 


523 


ESO 338-1G04 


267 (10) 


2.5e-12 


2.5e-12 


8.1e-13 


0.98 


1.69 


3.12 


15.8 


479 


SBS 0335-052 


3.40 (0.23) 


-4.0e-13 


3.6e-13 


9.9e-14 


-1.12 


-4.8 


3.62 


-27 


808 


NGC 6090 


147 (49) 


6.6e-13 


1.4e-12 


2.1e-13 


0.46 


82 


6.66 


62 


180 


IRAS 08339+6517 


157 (23) 


3.1e-12 


2.3e-13 


6.9e-13 


1.36 


3 


3.26 


45.6 


140 


Tololo 65 


22.7 (0.75) 


5.e-14 


1.8e-13 


4.9e-14 


0.28 


1.15 


3.6 


9.1 


1153 



Table 4. Integrated fluxes and equivalent widths for the six galaxies of the sample. The integration aperture is defined by masking 
regions below a certain threshold based upon H/3 flux (/h^j > 5cr, where cr is the background standard deviation in the same line). 
The aperture size is given in the first column in (")^ and in kpc^ in parentheses. The quantities are corrected for galactic extinction 



( Schlegel et al. 1998 ) but not for internal reddening, except for {Lya/Ha)c, which has been dereddened with Es-v.gda using Cardelli 
[et al.| ( |1989| ) parameterization. 



amount locally departs from the average. Including a large part 
of the Lya diffuse emission in the integration aperture could, for 
the same reason, contribute to this effect. 

We can retain from the above analysis that assuming simple 
dust extinction fails to recover the intrinsic LyafHa ratio, where 
the role of the dust is, in most cases, underestimated because of 
the resonant scattering phenomenon of Lye . 

4.2. Galaxy sample and Lya emission morphology 

We present here a comparison of our sample properties to 
their high redshift counterparts, such as Lyman Break Galaxies 
(LBGs). The FUV luminosity is that derived from integration 
within the apertures based on FUV background mask and the 
radius corresponds to Rav = yjareajTi. The galaxies span a 
large range in FUV luminosity 8.3 < logiLpuY / Lq) < 10.3 
that reaches, for two galaxies, the characteristic LBGs lumi- 
nosity, and are relatively compact systems, similar to LBGs 
radii range (log(R[/v(A:pc))~ 0-0.5). Accordingly, the FUV sur- 
face brightness Ifuv of the sample corresponds to that ob- 
se rved in Ultraviolet Lu minous Galaxies (UVLGs), as defined 

> 10^ Lq kpc~^ that clas- 
illustrates the 



by Heckman et al. (2005), with / 



sifies them in the "compact" category. Figure 10 



compacity of our galaxies among the UVLGs and LBGS, with 
their respective classification criteria on Lpuv and /puv overplot- 
ted. It is instructive to note, in reference to the present discus- 
sions related to local objects and the implications on high-z ob- 
servations discussed later, that Haro 11 and IRAS 08339+6517 
, could, regarding also their SFR and metallicity, be considered 
as UVLGs and LBG analogs. 



One of the consequences of the decoupling of Lya from non- 
resonant radiation is the diffuse emission halo observed in all our 



sample. We present in Fig. 1 1 the fraction of this diffuse com- 



ponent with respect to the total Lya flux as a function of the 
equivalent radius from the brightest emission source which has 
been isolated in each galaxy to produce the scatter plots of Sect. 
[3] The radius is yet derived from the surface size covered by the 
emission. SBS 335-052 and Tololo 65 are not represented on the 
figure since they do not show any bright (direct) emission that 
might be identified as photons production sources and are only 
dominated by a weak diffuse emission. At a distance of 2 kpc, 
the Lya emission reaches 40% of the total galaxy emission in 
only one galaxy, where for the remaining ones, it is still around 
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Fig. 10. FUV luminosity derived from F140LP flux (1500 A) in 
units of Lq versus galaxy radius (in kpc). The masking tech- 
nique based on the FUV background is applied to derive inte- 
grated Lpuv and the equivalent radius Ruv - y/area/n. The up- 
per dotted line corresponds to compact-large UVLGs threshold 
and lower ones are LBG limits in terms of Ifuv (see text). LBGS 
and UVLGs spaces are marked on the figure. 



20% or below. We need to integrate until 3 to 10 kpc (depending 
on galaxies) to account for the whole Lya flux. These consider- 
ations are relative to the threshold level reached by the observa- 
tions. We consider, because of the masking technique, that above 
80% of the total emission we are background limited. However, 
we can infer from these plots that the bulk of Lya emission con- 
sists of a low surface brightness region. 

The situation is slightly different in the galaxies not plotted 
for which the decrease of the surface brightness with radius is 
less steeper for reasons stated above. This illustrates the domi- 
nation of the photon diffusion mechanism in the escape of Lya 
photons, which appears at low surface brightness and extends to 
large physical scale (several kpc). 



4.3. Age and evolutionary effects 

We have computed the mean age for each galaxy from the SED 
fitting output in each resolution element with two stellar compo- 
nents as free parameters assuming in both of them an instanta- 
neous burst, and averaging over the whole integration aperture 
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Fig. 11. Lya emission morphology. The left plot shows the Lya 
surface brightness integrated in radius bins against the distance 
from the brightest emission center On the right plot, we see 
the evolution of the Lya emission fraction cumulated in varying 
apertures, as a function of the equivalent radius of these aper- 
tures. It shows what fraction of Lya is emitted within a given 
aperture around the bright emission center. SBS 335-052 and 
Tololo 65, which show no direct bright but only weak diffuse 
emission, are not represented. 

by weighting the age in each resolution element by the corre- 
sponding Ha luminosity: 



mean age 



Zi Lna, i X agej 



(2) 



except for Tololo 65, for which no HST Ha image is available 
and the luminosity in the B band is used for the weighting. We 
also carried out Monte Carlo simulations in order to estimate the 
errors on the weighted age. Each pixel has been resampled with 
1000 data points and the fitting procedure is applied to this new 
sample. The standard deviation gives the 1-cr errors plotted in 
Fig. [121 

The equivalent widths of strong hydrogen recombination 
lines are known to be, in principle, good age indicators since they 
measure the ratio between young ionising over old non-ionising 



radiations (Leitherer 2005 and references therein). In Fig. 12 
we observe an anticorrelation between Ha equivalent width and 
the age of the galaxy. This is what is expected by SED mod- 
els (starburst99 for instance) for both Ha and Lya equivalent 
widths, for ages > 1 Myr. This is due to the decrease of the ion- 
izing photons quantity and the increase of the number of stars 
that contributes to increase the continuum with time. However 
the dispersion observed is symptomatic of the complexity of this 
indicator in practise. EWua can be affected, among other effects, 
by the difference in reddening between nebular and continuum 
radiations and also the contribution of underlying older stellar 
population diluting the continuum. 

Plotting the Lya equivalent width yields a different result. 
EWi,ya does not follow this evolutionary sequence where an in- 
crease is observed as function of the age. The present obser- 
vations show the additional complexity of Lya with respect to 
Ha radiation, since we are comparing in the same plot different 
objects with different mechanisms at work, such as expanding 
shell in some ones and static medium in others, and/or probably 
clumpy ISM, that makes Lya quite difficult to interpret in a gen- 
eral way. Large variations in the observed EWi^ya are found in 
LBGs where most of them have likely the same intrinsic EWiya- 
[Schaerer & Verhamme] ( |2008| l have shown that for an extinc- 
tion of E(B-V)~ 0.3 an intrinsic emission with EWi^ya > 60 
A is transformed into an absorption, hence a negative observed 
EWiya- Radiation transfer and dust effects can lead to large dif- 
ferences between intrinsic and observed Lya equivalent widths. 



In addition, the apparition and the evolution of superwinds is 
yet function of the age, and we can expect that EWly^ rises in 
the presence of ISM kinematical effects. Because of the limited 
statistics, this could be considered only as a first step toward a 
more detailed and significant investigation. It is however worth 
noting that the small age (~ 2.6 Myr) of SBS 335-052 agrees 
with our discussion about this galaxy as a young starburst em- 
bedded in a static Hi leading to a damped absorption. This is also 
true for Tololo 65 (~ 3.3 Myr) where no direct emission is seen 
and for which GHRS spectrum shows a damped absorption. 
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Fig. 12. Age correlations: The left plot shows the integrated 
EWiya against the mean age integrated over the galaxy. On the 
right, EWua is plotted as a function of the age. 1-cr errors on 
the age fit are also shown. See text for age calculation and errors 
estimate. Each point represents a galaxy and the legend is given 
in the inset. 



4.4. Reddening correction and star formation rate 

The evolution of the Lfir/^^fuv ratio provides an alternative to 
evaluate the reddening estimate reliability. We show in Fig. 13 
a weak correlation between this luminosity ratio and the UV 
continuum slope p derived from the SED fit. We have plotted 
on the same figure (in blue) the predicted relationship (Kong 
let al. |2004) between Lfir/Lfuv and /3 following Meurer et al. 
( 19 99|l. dm galaxie s have the same dispersion behavior as those 

overplotted as red points, which are 
luminous in the IR. NGC 6090 falls 



of iBurgarella et al. (2005 



though more dusty and 
among the red points as it is classified as LIRG and is a very 
dusty starburst with -Eg-v^gas ~ 0.75). Three of our galaxies (NGC 
6090, Haro 1 1 and IRAS°08339-h6517 ) appear to be Luminous 
Infrared Galaxies (LIRGs, log(LFiR) > 1 1 Lq). In contrast with 
lUE starburst sample of Meurer et al. ( 1999 1, Xjoldader et al. 



(2002) found, as for our galaxies, that LIRGs and ULIRGs of 
their sample fall above the line. |Seibert et al.| ( |2005) have also 
derived, using a large sample of galaxy types, a correction that 
lowers the empirical Lfir/Lfuv-^S reddening relation and is thus 
also in disagreement with our observations. The observed dis- 
crepancy suggests that a simple empirical law, is actually not 
representative of the observations and the galaxies suffer appar- 
ently a higher attenuation than suggested by /5. Star-formation 
rate is an essential diagnostic tool of the evolution of galaxies. 
It allows, since the high-redshift universe has became reachable, 
to study star-formation episodes and their evolution on a wide 
range of epochs (e.g. |Madau et al.|1996[ Giavahsco et al.|2004[ l. 
Many indicators that ranges from rest-frame UV to infrared are 
used to estimate the SFR, and a commonly used one in the dis- 
tant universe is based on Lya emission, since the window opened 
on the high-z candidates ended up with using the Lya line as the 
main detection and investigation tool. One of the most critical is- 
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Target 




logiLFIR) 


SFR(UV) 


SFRCLyQ-) 


SFR(FIR) 


SFRCHff) 


(SFRt,v/SFRi,„) 


(SFRu v/SFRlto )cor 


Haro 1 1 


10.3 


11.1 


5.1 


1.07 


21.7 


16.33 


4.77 


0.53 


ESO 338-IG04 


9.76 


9.8 


1.56 


0.3 


1.1 


3.9 


5.20 


3.20 


SBS 0335-052 


9.1 


9.4 


0.34 




0.4 


1.1 






NGC 6090 


10.03 


11.4 


2.9 


1.1 


43.3 


20.7 


2.64 


0.015 


IRAS 08339+6517 


10.3 


11.0 


4.4 


2.2 


17.3 


14 


2 


0.85 


Tololo 65 


8.3 


8.4 


0.05 


0.008 


0.05 


0.24 


6.25 


2.6 



d525A. li is the flux density in F140LP filter. 



Table 5. Star formation rates and luminosities. Lfuv is calculated with /Ix/^i where A- 

LfiR is from Ostlin et al. (2008 1. Both are expressed in units of Lq^boI (S.SxlO^'' ergs s"'). Star formation rates are derived from 
the integrated fluxes over apertures based on sky background threshold and using the calibration of Kennicutt (1998) (presented in 
units of Mq yr"^). All quantities are corrected for galactic foreground extinction. The last column ratio only is corrected for internal 
reddening using ZiB-Kstais for SFRyv correction and EB-v.gas for SFRLya one. 




Fig. 13. L/R to Lfifv ratio versus the UV slope j6. Our galax- 
ies are plotted with the same symbols as in previous figures. /3 
is determined from the SED fitting. The remaining red points 
are galaxies from ELAIS SI observations in Burg arella et al.| 
( [2005 | l. The law describing the relationship between FIR to FUV 
luminosity ratio and the slope /3, deduced from UV observations 
( [Kong et al.|2004| l is overplotted (dashed blue line). 



sues related to estimate the SFR, from UV or optical indicators, 
is the correction for internal reddening. One needs to properly 
estimate the dust obscuration and therefore the intrinsic flux to 
make a correct conversion to star-formation rate. This issue is 
even more critical for the indicator based on Lya emission, con- 
sidering the radiative transfer complexity of this line discussed in 
this paper Table |5] summarizes the SFRs computed from differ- 
ent indicators using |Kennicutt| ( |l998| l calibrations. The conver- 
sion from flux to SFR from these calibrations assumes a contin- 
uous star formation regime in the equilibrium phase, whereas in 
our SED fitting procedure, used for the Lya continuum subtrac- 
tion or age estimation, we have assumed an instantaneous burst 
scenario. Nevertheless, translating our results into star-formation 
rates gives a useful comparison with previous works using these 
widely used calibrations. 

UV emission from galaxies traces the young stellar popula- 
tion and a conversion from UV luminosity to SFR can be com- 
puted across the UV range (1250-2500 A). This estimation is 
very sensitive to dust attenu ation, beca use of the wavelength do- 
main and the patchy ISM ( ,Kennicutt 1998| '] Alternatively, nebu- 



lar emission lines, such as Lyor, provides information on the ion- 
ising flux of the young massive stellar population with lifetimes 
< 20 Myr Therefore, it provides a quasi-instantaneous estima- 
tion of the current star-formation rate. For case B recombination 
theory ( |BrockIehurst| 1971] ) we have: 



SFRLy„(Mo yr-') 



9.1 X lO""*^ L(Lya) (ergs s"^) 



(3) 



Yet, in addition to the IMF dependency, this method is highly 
sensitive to the extinction correction as we have mentioned be- 
fore. F igure [l4] is illustrative of the reddening correction issue. 
Figure 14(a) represents SFR(LyQ') versus SFR(UV): for dark 



points neither SFR is corrected and for red points SFR(LyQ') 
is corrected using EB-v.gas and SFR(UV) using ^B-v^stais ■ In 
the purpose of comparison, data from literature are also in- 
cluded. It is interesting to note the consistent discrepancy be- 
tween SFR(LyQ') and SFR(UV) (dark points) for both our sam- 
ple and objects from literature, including low- (Giavalisco et al. 
19961 and high-redshift galaxies ( [Taniguchi et al. 2005 Ajiki 
et al. |2003[ [Tapken et al. |2007] l. The lack of points under the 



line of equality at low SFR(LyQ') for high-z observations, is a 
consequence of the completeness limit, for LAEs in particular, 
excluding faint Lya emitters. Indeed, we expect, because of the 
resonance effects of Lyo- to observe a scattered distribution be- 
low the line of slope unity providing one uses a statistically 
significant sample. Such a distribution would help to charac- 
terize a potential upper limit of the discrepancy between res- 
onant and non-resonant SFR indicators. Such a diagram could 
also serve as an important probe for the galaxy evolution from 
the damped and/or young systems lying at SFRLyt, ~ 0, to more 
evolved starburst events with higher ionised gas fraction and/or 
undergoing feedback outflows hence approaching the line where 
SFRLyff ~ SFRuv- 

For our galaxies, in the case of observed values, SFR(Lya') is 
systematically below the equal values line, underestimating SFR 
by a factor between 2 to 6 with respect to SFR(UV). This dis- 
crepancy, usually observed in high-z galaxies also, emphasizes 
the highest attenuation of Lya emission line with respect to UV 
continuum. The correction for dust attenuation, one of the main 
responsible of this discrepancy, brings the two star-formation 
rates in a better agreement except for NGC 6090. 

When SFR derived from the UV is not corrected for absorp- 
tion, the SFR measured from the infrared should, in principle, 
be complementary, since radiation that is stro ngly ab sorbed in 



shows the 



the UV is re-emitted in the thermal IR[^ Figure 14(b) 
"total" SFR (SFRuv + SFRfir) where SFRf/y is not corrected. 



SFRuv(Mo yr-') = 1.4 x IQ-^** L,(ergs s"' Hz-') 



SFRfir (Mo yr"') = 4.5 x 10-^* L™(ergs s"') 
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Fig. 14. Star formation rates: Figure 14(a) shows SFR derived from nebular emission line Lya versus SFR derived from UV contin- 
uum. The dark points represent undereddened SFRs. Note that because of the logarithmic scale and for the sake of readability, we 



have put galaxies for which SFRj,yQ,-0 at SFR^yg- O.OOl. Observed SFRs from l iterature at low-z ([G iavalisco et al.|1996 and lUE 
- •™ '■ ~ '-'^^^^^ — n^TTTT:^, „ . . „ , J— = ,yi ]•■ — ■ 



data from NED) and high-z ( [Taniguchi et al.|2005[|Ajiki et al.|2003 ], for LAEs) ( jTapken et al.|2007[ for LBGs) are also overplotted 



and explained in the legend. For our sample, the red points represent SFRiya dereddened using Eb^v.^rh and SFRyv dereddened 
using EB-v,stais- The dashed line is for SFR^,^ - SFRj/y. Figure 14(b) represents SFR(UV-hlR) (where SFRj/y is not corrected for 



extinction) versus different corrected SFRs. The legend is also marked on the figure. Units are in Mq yr 



versus dereddened SFR^,,„ and SFR/za with Es_v,gds (dark and 
blue points), and SFRuv with Es-v.stias (red points). The cor- 
rected SFR(LyQ') is still below the total SFR for most of the 
galaxies except for NGC 6090 where the dust correction highly 
overestimates the total SFR (we have discussed possible reasons 
for this overestimation in Section 4.1 1. Conversely, SFR derived 



from corrected Ha (blue points) luminosity yields a better result. 
Similarly, the dereddened UV estimator (red points) places the 
galaxies rather close to the line of equality. 

Finally the different results found when deriving SFR from 
Lyo- and others indicators such as Ho- or UV, and in particular the 
failure of Lya indicator to recover the total SFR (UV + IR) even 
when corrected for reddening, are indicative of the decoupling 
of resonant Lya and non resonant (e.g. UV continuum or Balmer 
lines) radiation with respect to the dust obscuration and therefore 
the difficulty to use this line as a reliable star formation indicator 

As an alternative, we can take advantage of the available in- 
formations of our observations to make a better estimate of the 
star formation rate when only Lya luminosity is known. We can 
calculate the Lya escape fraction using the corrected Ha flux 
and assuming the case B recombination ratio Lya/Ha: 



Uc(Lya) = f(lya)/(^Jxf(Ha)c) 
f{Ha)c = f(Ha) x io(1"48x£(b-v).,,) 



(4) 
(5) 



where f(lya) is the observed flux and /(Ha)c is the Ha flux cor- 
rected pixel by pixel for internal reddening using |Cardelh et al. 
( 1989 1 extinction law. Unlike high-z observations, where in the 
best case only the global UV slope is known, our present study 
gives access to accurate extinction information. Therefore, fg^c is 
a good estimate of the intrinsic Lya flux, since it takes into ac- 
count both dust obscuration and resonant scattering mechanism 
and hence allows us to correct the SFR/,yQ, calibration (Eq. |3|l as 
foUows: 

SFRLy,(Mo yr-') = (l//,,,)x9.1 x 10-^^ L(Lya) (ergs s"') (6) 



The escape fractions obtained for the six galaxies are reported in 
Table |6] We can assume, in general terms, an escape fraction of 
5% as a mean "statistical" value to deduce the appropriate SFR 
when fgsc is not available. This correction is evidently subject 
to uncertainties due to f^sc variations, but is in any case more 



representative of the reality than standard calibrations ( Kennicutt 
|1998| for instance). 

The revised star formation rate based on L(Lya) is equiva- 
lent, according to the definition of the escape fraction in Eq. |4] 
to SFR(Ha) corrected for reddening. Hence, it is represented on 
the right plot of Fig. 14 by the blue points, which appears to give 
a better estimation of the total SFR than that given by the SFR/,yQ. 
corrected for only dust obscuration. 

4.5. Implications for high-redshift galaxies 

We show in figure[T5]how fgsc is function of the extinction. We 
observe that the Lya escape fraction is somewhat correlated with 
^B-v.gas but with an important dispersion in, moreover, a quite 
limited sample. At a galaxy scale we see the manifestation of 
the dust obscuration effect on emergent Lya radiation, although 
in a purely dust regulated model, this correlation would have 
been more striking. It is also interesting to note that correcting 
the total Ha flux using the mean Es-v.gas lead to higher escape 
fractions by 15 to 40 % that those calculated using pixel level 
corrections. In high redshift observations only global correction 
is possible in absence of spatial resolution. 

The rate of escaping Lya photons does not exceed 10% in 
this sample. In addition, we have outlined the ubiquitous halo 
of diffuse emission present in all the observed galaxies, and we 
are able to get a feeling of how important is this contribution 
in section 4.2 The diffuse component represents the bulk of the 
Lya emission (Table |6]l. It is very likely that such low surface 
brightness emission remains undetectable in high-redshift galax- 
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ies. This point out the caution that should be taken when de- 
riving physical quantities, such as star formation rates (SFRs), 
from Lya alone, since we are dealing with only a small fraction 
of escaping photons, and we are probably missing the major- 
ity of this fraction. Simple extinction correction of SFR could 
also be at fault since such method has failed to recover the value 
determined by the recombination theory in most cases and the 
behavior of resonant radiation according to the dust content is 
unpredictable without any complementary information (HI dis- 
tribution, gas kinematics, etc.). However a better, and more re- 
alistic estimate could be obtained using the calibration proposed 
in Eq. |6] 



prompt a further investigation with an extended sample of galax- 
ies and complementary spectroscopic study. 
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Fig. 15. Escape fraction of Lya photons (in percent, see text 
for details on the f^^c determination) as a function of extinc- 
tion E(B-V) in the gas phase. We have included the two compo- 
nents (knots A and B) of the interacting system NGC 6090 and 
treated them separately by a masking procedure. We observe a 
decline in the amount of escaping photons when increasing the 
dust amount. We note also that most of the galaxies have a small 
fesc around 3% or below. 



Because of the difficulties discussed above, high-z star for- 
mation rates based on Lya are generally underestimated with re- 
spect to that derived from UV for instance. Discrepancies have 
been observed between the two estimation methods where SFR 
based on the Lya luminosity was smaller by a factor of two or 
more than that based on the UV continuum ( |Hu et aL]|2002[ 
Kodaira et al.|2003| l. [Tim"guchi et al.| ( [2005) found that SFRs de- 
rived from Lya for their sample of z ~ 6.6 LAEs lie a factor of 5, 
on average, below those based on UV continuum. T apken et al.| 
( |2007| l found similar discrepancy for their UV-selected galaxies. 
This is also what is seen at low redshift where a factor of 2 to 6 



is found (the present work. Section 4.4 1. 

Lya emission line has became a powerful tracer of star for- 
mation at high-redshift. However, it is clear that using only the 
Lya luminosity or equivalent width for the characterization of 
star formation episodes leads to many uncertainties which inhere 
in the Lya radiative transfer complexity and the uncertainties of 
the SFR calibration methods. Despite differences due to obvious 
evolution effects between local and high-z star-forming galaxies 
(see Mas-Hesse et al. 2003[ ), the physical processes governing 
the radiative transfer and the escape mechanism of Lya pho- 
tons should be the same, and motivate the present extrapolation. 
Nevertheless, we stress again the limited number of galaxies in 
our sample that precludes any statistically significant study and 



Target 


f 


^diffuse 


E(B-V) 


E(B-V) 


AGE 




(%) 


(%) 


gas 


stars 


(Myr) 


Haro 1 1 


2.6 


74 


0.33 


0.07 


11 


ESO 338-IG04 


8 


70 


0.08 


0.08 


5.4 


SBS 0335-052 


0.25 


~ 100 


0.21 


0.04 


2.6 


NGC 6090 


0.56 


73 


0.76 


0.18 


25 


IRAS 08339+6517 


7 


65 


0.12 


0.03 


18 


Tololo 65 


1.7 


~ 100 


0.2 


0.12 


3.3 



Table 6. Integrated properties and Lya emission characteristics, 
quantities are integrated over the same apertures defined in table 
|4] The first column contains the escape fraction of Lya photons 
for each galaxy assuming a case B Lya/Ha recombination ra- 
tio (see text for details). For SBS 0335-052, only positive Lya 
contribution is used to derive an upper limit for f^sc- Second col- 
umn represents the contribution of the diffuse component in the 
whole Lya emission determined by masking the bright "direct" 
emission sources if any. The mean nebular extinction fig-v^gas 
is derived from the ratio of integrated Ha to H/? fluxes and is 
showed in column 3. The fourth gives the continuum extinction 
determined from the SED fitting procedure. Last column shows 
the age issued from the SED fitting and averaged on the total 
galaxy pixels following Eq. |2] 



5. Conclusions 

Combining space (HST) and ground-based (NOT and NTT) 
observations we have mapped the Lya emission and the dust 
content in six nearby star-forming galaxies. We have com- 
pared the extinction E(B-V) produced from the Balmer decre- 
ment Ha/H/3 to several parameters like Lya emission, equiva- 
lent width or recombination ratio Lya/Ha at a small scale in or- 
der to disentangle the role of the dust from other parameters. 
Implications for high-z studies, via the global properties of the 
galaxies have also been investigated. 

Our galaxies exhibit different Lya morphologies from emis- 
sion to damped absorption or combination thereof: 

- In systems with emission and absorption (namely Haro 11 
and ESO 338-04) we found Lya photons emerging from re- 
gions with similar or even higher extinction than those where 
Lya is seen in absorption. We point out the role of the ISM 
distribution, where in the case of clumpyness morphology 
Lya photons escape preferentially to Ha ones leading to an 
observed Lya/Ha ratio higher than the theoretical level cor- 
rected for extinction. 

- In objects that show no strong absorption (NGC 6090 and 
IRAS 08339+6517 ) we observe no clear correlation be- 
tween Lya and the dust content. The Hi kinematics may play 
a more significant role in the escape of Lya photons as con- 
firmed by kinematics studies which have shown large ISM 
outflows in both systems. 

- SBS 335-052 is a Lya absorber with a great Hi column den- 
sity coverage which is believed to be static with respect to 
the emitting region. We estimate an age (< 5 Myr) in agree- 
ment with the picture where the starburst is too young to have 
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ionised the surrounding gas or driven an outflow. We observe 
precisely what is expected from the resonant nature of Lyor 
in a static neutral gas: a damped absorption with a declin- 
ing relationship between Lya and E(B-V) indicating that the 
dust is, in this case, the main regulator of Lya escape. 

When investigating global parameters of our sample, we 
found that simple dust extinction correction fails to recover the 
intrinsic Ly afHa ratio, where the role of the dust is, in some 
cases, underestimated because of the resonant scattering, and in 
other cases, overestimated because of the clumpyness distribu- 
tion of the ISM. We observe neither no evident correlation be- 
tween EWijya and the reddening. The observed Lyor escape frac- 
tion is found to not exceed 10% in our sample and is, for most 
of our galaxies, about 3% or less. 

The resonant decoupling of Lya from non-resonant radiation 
leads also to an ubiquitous diffuse halo with low surface bright- 
ness. It represents the bulk of the Lya emission and extend to 
regions at several kpc from emitting regions and which are not 
reached by Ha or continuum radiation yielding high EW^ya- 

Because of the radiative transfer complexity of the Lya line, 
star formation rate (SFR) based on Lya lead to different results 
than SFR derived from other indicators (from UV for instance), 
and fails to recover the total SFR (UV + IR), and this even 
when corrected for dust obscuration, preventing any determina- 
tion of the intrinsic star formation rate. We propose therefore 
a more realistic calibration for the SFR when information on 
Lya only is available (which is usually the case for high-redshift 
surveys), accounting for dust attenuation and resonant scattering 
phenomenon through the Lya escape fraction. 
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